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In this paper we compute the coefficients of the heat kernel asymptotic expansion for 
Laplace operators acting on scalar functions defined on the so called spherical suspension (or 
Riemann cap) subjected to Dirichlet boundary conditions. By utilizing a contour integral 
representation of the spectral zeta function for the Laplacian on the spherical suspension we 
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find its analytic continuation in the complex plane and its associated meromorphic structure. 
Thanks to the well known relation between the zeta function and the heat kernel obtainable 
via Mellin transform we compute the coefficients of the asymptotic expansion in arbitrary 
dimensions. The particular case of a d-dimcnsional sphere as the base manifold is studied as 
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well and the first few heat kernel coefficients are given explicitly. 

I. INTRODUCTION 

The study of the heat kernel for Laplace type operators has been a topic of major interest in 

> 

the literature. The reason lies in the fact that it represents a very versatile tool for the analysis of 

\^ . various problems. Its interest in mathematics resides, especially, in geometric analysis thanks to 



the well known relation between the heat equation and the Atiyah-Singer index theorem 
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In physics, the heat kernel of an elliptic operator and its associated spectral zeta function are 
particularly useful in the ambit of quantum field theory in order to compute the one-loop effective 



action and the Casimir energy 
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The computation of the heat kernel, however, is not a trivial task. In fact, its exact form can only 
be found for highly symmetric backgrounds when the spectrum of the operator is explicitly known. 
For this reason, various approximations have been developed and used throughout the years. One 
of the most important approximation schemes is the short time asymptotic expansion of the heat 
kernel also known as Schwinger-DeWitt expansion 15l4l8l|. Very well established methods have 



been developed in order to compute the coefficients of the heat kernel asymptotic expansion for 



elliptic operators on compact manifolds with and without boundary (see e.g. [a. I la. |29i. l30l. l35l. |4(| ) . 



A very efficient way for evaluating the coefficients on a manifold with boundary relies on the use 
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of the spectral zeta function for special cases. This is done by exploiting the wel 

m 



known relation 



351 ] . We provide 



between the trace of the heat kernel and the associated spectral zeta function 
a new example here that can be used for this purpose. 

It is well known that for second order self- adjoint elliptic partial differential operators on com- 
pact smooth manifolds the spectrum X n is discrete, bounded from below and it forms a mono- 
tonically increasing sequence of numbers tending to infinity. The spectral zeta function for such 
operators is defined as the following sum 

oo 

c(*) = E A « s ' ( L1 ) 

n=l 

which is convergent for 3ft(s) > D/2, with D being the dimension of the manifold under considera- 
tion. One can analytically continue C( s )> i n a unique way, to a meromorphic function in the whole 
complex plane which coincides with (jl.ip in its domain of convergence. In this paper we will utilize 
spectral zeta function techniques in order to evaluate the coefficients of the heat kernel asymptotic 
expansion for a Laplace operator, with Dirichlet boundary conditions, on a singular Riemannian 
manifold known as the spherical suspension. A related work on the spherical suspension is [241 ]. 
where the functional determinant for a Laplace operator and the conformal anomaly have been 
analyzed. 

The current paper represents somewhat a generalization of the conical geometries studied in 
[7]. The difference is not in the mathematical technique used, but in the radial eigenfunctions of 
the Laplace operator that one finds. In the case of the generalized cone, a prototypical singular 
Riemannian manifold, and in many other ones considered in the literature, the radial eigenfunctions 
are Bessel functions. Well developed methods have been described in order to use their uniform 
asymptotic expansion for the analytic continuation of the spectral zeta function. In the case of the 
spherical suspension considered here the radial eigenfunctions are associated Legendre functions. 
To our knowledge, the analytic continuation of the associated spectral zeta function by making use 
of the uniform asymptotic expansion of the Legendre function is analyzed here for the first time. As 
we will see in the next sections, the technical details are quite different from the Bessel case which 
makes this work particularly interesting. It is worth mentioning that zeta function regularization 
techniques involving associated Legendre functions have been considered in the study of the one- 
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loop effective action for gravitons in a DeSitter spacetime [5|. In that paper, however, a different 
procedure was used from the one presented in our work. 

The layout of the paper is as follows. In the next section we will describe the geometry of the 
spherical suspension and present the eigenvalue problem we are interested in. In the subsequent 



sections we will construct the zeta function in terms of a contour integral and explicitly find its 
analytic continuation. This will allow us to compute, in section |Vj the heat kernel coefficients on 
the spherical suspension in terms of those on its base. As a particular case in section EH we study 
the situation in which the base manifold is a <i-dimensional sphere which will permit us to find 
more explicit results. 

II. THE GEOMETRY AND HEAT KERNEL 

The manifold that we will consider in this paper is termed spherical suspension or Riemann 
cap. This is defined as the D = (d + l)-dimensional manifold ^i = I x ,JV where X C [0,#o]> 
with #o S (0, 7r), and where JV represents a smooth, compact Riemannian ci-dimensional manifold, 
possibly with boundary, which will be referred to as the base manifold. The local geometry of ^# 
can be written in hyper spherical coordinates as [7j, |24| 

ds 2 = d9 2 + sin 2 6dY? , (2.1) 

where dT, 2 represents the line element on the base manifold jV . It is easily seen, from (|2.ip . that 
the Riemann cap ^ in general possesses a singularity and, therefore, it belongs to the class of 
singular Riemannian manifolds. Since sin#o ~ #o i n the limit of small angles 9q the metric (|2.1|) 
reduces to the one representing a generalized cone [7]. This special limiting case will provide a 
robust check of the results that we will obtain in this paper. 

In hyperspherical coordinates the Laplacian on ^, namely A^, can be written as follows 

d 2 d 1 

A ^ = W + dcote de + ^o A ^ M 

where A^ is identified with the Laplacian on the base manifold JY ' . The eigenvalue equation 

- A.^</? = a 2 tp , (2.3) 

with if belonging to the space £ 2 (^#), is separable in this system of coordinates and a general 
solution to (|2.3p can be written as a product 

( p= (sin 6)^ ?P(6)H . (2.4) 

Here, the functions % are defined to be the harmonics on ,jV having degeneracy d(\) and satisfying 
the following eigenvalue equation 

- \_, Y U = \ 2 U . (2.5) 



As will become clear later, the eigenvalue problem resulting from 



d , 
-A^r + -j J <P = u V, 



(2.6) 



has certain advantages compared to (|2.3p . In this case the function ^ in (j2.4|) obeys the differential 
equation of Legendre type 
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(2.7) 



I / sin* 

This is the eigenvalue problem we will concentrate on and we will make remarks at appropriated 
places on how the original eigenvalue problem f|2.3j) can be recovered and what the additional 
complications are. 

Let us next turn our attention to the associated heat kernel and its asymptotic expansion. For 



explanatory purposes, let M be a smooth manifold with boundary dM. It is well 
£ 2 -trace of the heat kernel has the following small t asymptotic expansion 



mown that the 
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where the coefficients 



Ek 
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dV et + dS fk , 

M 2 JdM 



(2.8) 



(2.9) 



consist of a volume and a boundary part, with ek and fk being polynomials in the local curvature 

2 2 

of the manifold. There exists a deep relation between the trace of the heat kernel and the associated 
spectral zeta function, through the Mellin transform, which is provided by the formulas 
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with s = n/2, (n - l)/2, • 
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= r(s)ResC(s) , 

(2Z + l)/2 for/ GN , 
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-C(-p) , 



(2.10) 
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where p £ No. These relations produce any coefficient of the heat kernel asymptotic expansion 
in terms of either the residue or the value of the associated zeta function at a specific point and 
they will be used in this paper for the computation of the heat kernel coefficients on the spherical 
suspension. It is clear that the heat kernels associated with the eigenvalue problems (J2.3J) and (12.6|) 
are related by a factor of e~ d 4 ' 4 . So, as far as the heat kernel is concerned, considering (j2.6|) is as 
informative as considering (12. 3D . The corresponding changes in the heat kernel coefficients can be 
trivially obtained. However, for a full analysis of the associated zeta functions the presence of a 
mass term gP/4 would cause major complications and we will further comment on these later. 



III. THE SPECTRAL (-FUNCTION 

Let us now begin to study the spectral zeta function of the eigenvalue problem (12 .6p . By 
introducing 



M= \/A 2 + ^t^> (3-1) 



a solution <j> to the differential equation (|2.6p which is regular for 9 — > is 24] 



^(0)=p-? (cosfl), (3.2) 
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where P_w 2 _i_ w corresponds to the Ferrers representation of the Legendre function 
The spectral zeta function for the problem under consideration is defined as follows 

C( S ) = J>- 2s , (3.3) 

where we will assume that no negative eigenvalues occur so that we can use the standard branch cut 
of the logarithm. Since jV for now is unspecified, we will express the spectral zeta function on the 
whole manifold in terms of the zeta function C.yy{ s ) associated with the operator — A t/ y + (d— 1) 2 /4. 
Its definition is [L3J 

0(s) = J>(/^- 2s . (3.4) 

Although the base manifold has not been explicitly chosen, one is still able to impose boundary 
conditions. For definiteness we will impose Dirichlet boundary conditions on the solution (|3.2p at 
9 = 9$. By doing so, we are led to the relation 

p :i/ 2+ > os ^) = o. (3.5) 

The last equation determines implicitly the eigenvalues lu 2 . 
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It is important to make a remark at this point. It is well known [HJ, I 111 . 
problems, such as the one considered in this paper, the heat kernel asymptotic expansion possesses 
a non-standard logarithmic term. This means that instead of a small t asymptotic expansion of 
the form (12.81) one has 



JC(t) ~ J2 Agrr- + & In t . (3.6) 



k=0 

This unusual behavior of the heat kernel leads to the appearance of a pole in the spectral zeta 
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function of ^# at s = and, in addition, one can prove that <$ = — Res£(0) [7J]. It has been shown 
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2J j that for the spherical suspension the following relation holds 
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-ResC(O) , 



(3.7) 



which is the same that one obtains in the framework of the generalized cone [7]. It is clear that 
the coefficient of the logarithmic term in the heat kernel asymptotic expansion can be explicitly 
evaluated once the residue of C/K s ) at s = — 1/2 is known. 

In order to compute the spectral zeta function, which will be used to evaluate the heat kernel 
coefficients, we will temporarily consider the massive case. This means we will add a constant term 
m 2 to the eigenvalues uj 2 , which at a later stage will be sent to zero. The reason for proceeding 
in this way is twofold. First, the presence of the mass parameter will allow for an easier analytic 
continuation of relevant integrals. Second, we will obtain an expression for the spectral zeta function 
which is valid also for a massive scalar field, results useful when Casimir energies for massive fields 
in the given setting are considered. 

To proceed, we express C[s) in terms of a contour integral in the complex plane, valid for 



»(«) > (d+l)/2 as 
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C« = £ d(n) JL jf dk (k 2 + m 2 ) - s A In F-J 1/2+k (cos O ) , 



(3.8) 



where T represents a contour that encircles all the roots of ~P~i/2+k ( cos ^°) on ^ ne P os iti ve rea l ax is 
in the counterclockwise direction. By deforming the contour of integration T along the imaginary 
axis and by using the property P^ M = P_i_„ [22J, |37J], the expression (|3.8|) can be rewritten in the 
following form 
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Here, the integral is well defined in the strip 1/2 < 3?(s) < 1 and £(s) above has to be interpreted 
as the function resulting from analytic continuation out of that strip. By conveniently changing 
the integration variable z — > (iw/ sin #o , one can rewrite (|3.9p as 

s d 



((s) = Y,d(n) 
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H 2 w 2 
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(3.10) 



By following the technique developed in pHS], we now utilize the uniform asymptotic expansion 
of the Legendre functions for jjl — > oo. By exploiting the WKB method, one obtains the following 



uniform asymptotic expansion valid for /j, — > oo and fixed w 
P~l A _ A kw (cos O ) " 
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(3.11) 



where we have introduced the notation 



and 
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1 + ln 



IV 



VYT^' 



w 



COS 6q 

VTT^ 2 ' 



arctan 



smfo 
w 



.cos 9 + Vl + w 2 } sin »o 
The functions x J f ;(i / ) in expression (13. llf) are defined by the formula 
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(3.12) 
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where B21 are the Bernoulli numbers and the functions & n {v) satisfy the recurrence relation 24], [33] 
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with $o = l an d 7 = w/ svo-Oq. For the purpose of this paper we will actually need the uni- 
form asymptotic expansion for the logarithm of the Legendre function (|3.1ip . To this end, it is 
straightforward to obtain 



lnP_f „ w (cos 9 ) ~ In 



2 sin fQ 
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(3.16) 



with the newly introduced functions £l n { v ) being defined through the cumulant expansion 
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The first few Q. n (v) are listed for completeness in appendix El 

By adding and subtracting N leading terms of the asymptotic expansion (|3.16p , one can write 
the spectral zeta function ((s) in the following way (cf. [7H9l.l35l]) 

N 

C(s) = Z(s) + J2 Ms) , (3.18) 
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where the terms Z(s) and Aj(s) are defined as 
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(3.19) 
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(3.20) 



A (s) = Y / d(ix) 
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and finally 
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It is very important to notice that Z(s) is an analytic function in the region Sft(s) > (d— 1 — N)/2. 
This restriction is found by considering the behavior on the integrand in /U, which is fi~ 2s ~ and 
by recalling that the spectral zeta function on J/ converges for 3ft(s) > d/2 according to Weyl's 
estimate. This means that in any dimension D one can choose a suitable N such that Z(s) does 
not contribute to the residue of the zeta function £(s) in the region K(s) > (d — 1 — N)/2, and 
hence to the value of the heat kernel coefficients (cf. I2.1U|) [8|, 135 ]. 



IV. EVALUATION OF A_ x (s), A (s) AND A^s) 



In this section we are mainly interested in the non trivial task of explicitly evaluating the 
functions in (13.20h -( |3.22|) . According to the remarks made earlier, these are the only terms that 
will contribute to the coefficients of the heat kernel expansion. Let us start with the leading term 
j4_i(s). By recalling the expression for r] in (|3.13|) one obtains 
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(4.1) 



It does not seem possible to find a simple analytic expression for the above integral, however some 
progress can be made by considering the following expansion 
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The expansion above is obtained by first using the relation 
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(4.3) 



and then by expanding the arcsine for small arguments. By using the result (|4.2p in the integral 
(I4.1D . and by making the change of variable w — > ms'mOo/dit) one obtains 
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The integral in the above expression can be evaluated in terms of hypergeometric functions, and 
its explicit form is 
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which is well defined in the strip 1/2 — j < K(s) < 1. In order to be able to perform the sum over 
the eigenvalues \i we utilize the Mellin-Barnes representation of the hypergeometric function 31 ] 
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where the contour is chosen in such a way that the poles of T (u — k — ^) T (u + s + j — ^) he to 
the left and poles of r(— t) lie to the right of 'rf. 

In order to proceed, we shift the contour of integration ^ so that W€ < 0. During this procedure, 
the contour must cross the poles of T (u — k — |) on the real axis for u > 0, which are positioned at 
u = k + 1/2 with k £ Nq~ (as illustrated in Figure 1). By taking into account the above mentioned 
poles and the poles to the left of the contour, one obtains using the residue theorem 
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By substituting the last expression in formula (|4.5p . one can prove that potentially divergent terms, 
in m and sin#0) cancel exactly giving the final result 
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Figure 1: The figure illustrates how to deform the contour as done in passing from formula (|4.5j) to formula 



where in order to obtain the last formula we have used 



T(-x-n) = r _ n n+l TT _J_ 

r(-x + i) l ; l = {x + k 



(4.8) 



which is valid for x ^ —n and n G Nq", and can be proved by induction. 

The next term that we need to consider is Aq(s). The evaluation of the integral in (|3.2ip is 
rather standard (see e.g. [21, la]) and its explicit form is 
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The last term that we want to study is Ai(s) constructed from the functions Qi{v). For the 
computation of this integral one needs to find a general expression for Qi(y) in terms of w. By 
introducing the function, defined for i G N + , 

1 + (-!)* 
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where [x] represents the integer part of x, we can write the functions £li(v) as follows 
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The numerical coefficients Zq and zi &j can be obtained from the relations (J3.15P and f)3. lTj) 
and the polynomials Di(x) are the ones associated with the uniform asymptotic expansion of the 
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modified Bessel functions (see e 
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defined through the cumulant expansion 
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We would like to remind the reader that the functions Q{(v), up to the fifth order, are listed in the 
appendix lAl 

By taking into account the specific form of £li(v) we can write Ai(s), for convenience, as a sum 
of three terms 
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with the definitions 
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The first two integrals can be evaluated by following the techniques utilized in [7J,|a|. After a rather 
standard computation one obtains the results 
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The computation of the last integral, namely 7^3 (s), is actually more involved. Integrating by 
parts and performing the change of variable w — > msin#o/(^v^) one can recast (J4.18P in the form 

\r^ ,( \ sin 7rs \p y^ Zjjjj (cos 8 y +2b 
v J =1 b=x(i) 



lids) 



-sm 



-2s 



/'' 



,2 e ,; T ,2 fl _\-&-| 



Jo V M 2 



4 + 7? 



(4.21) 



which is valid in the region —j — b — i/2 < K(s) < 0. By exploitin g th e following general integral, 
valid for 3?(a) > and 3? (6) > 0, with non-vanishing real A and I? 22J, l3l| 



2\~<* 



(I 



dxx a - l {l-x) b - l (x + A 2 ) (x + B 2 ) 
A -2c B -2 d mm Fl L c,d,a + b;--L,-±), 



r(o + 6) 
where -Fi represents the hyper geometric function of two variables, we obtain the expression 

■>e y +2b ^ 2s T{s + b+j + t 



(4.22) 



IM 



E«E E 

M J=l 6=x(») 



: m 



r(s)r(6 + j + f) V^sin^o 



/' \ f£\ 2 ° 

m) 



x Fl{s + b + J + ^ b + ^ b + J + ^-^F ^^ 2 



(4.23) 



The hypergeometric function of two variables that appears in the previous equation can be ex- 
pressed in terms of a hypergeometric function of one variable thanks to the relation 



22 



m 



Fi(a, 6, b', b + b'; x, y) = (l- x)- b (l - y)- a+b 2 F 1 lb + b' - a,b,b + b'; 



y-x 

1 — X 



(4.24) 



The useful formula (j4.24[) allows us to write (I4.23P in the form 



ii,3(«) = -E d (^ 2s 'EE 2i,6j(cos0 o 



,i+26 



T(s + b + 



Lt&h 



3=1 6=x(0 



r( s )r(6 + j + i) 



+ 



m 2 sin 2 6>o^ " 



/< z 



x 1 + 



m 



2\ -J- s 



2 Fi -s, 6 + -, 6 + j + -; cos z 6<o ( 1 + 



2 ■ 2 /) \ ~l x 

m sm pq n 



(4.25) 



i 1 J \ A z V A 4 " 

The argument of the hypergeometric function in (J4.25P is smaller than one which guarantees that 
an expansion is applicable. As we have done before, we utilize a Merlin-Barnes contour integral 
representation for the hypergeometric functions appearing in (|4.25p . In this particular case, since 
the argument of the hypergeometric function is smaller than one, we close the Hankel contour to 
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the right. The resulting integral can be computed with the residue theorem where the poles of the 
integrand are located at the positive integer points. This yields an expansion in terms of positive 
powers of the parameter msin#o- On doing so, after some manipulations of the gamma functions 
and an expansion in powers of m of the term (l + m? sin 2 ^o/m 2 ) ' 3 (l + m 2 /^ 2 ) , we obtain 
the final result for 1^3 (s) as follows 

i i3 ( s ) = - J_ y y Zib {cos e y^ b r ^ s + b ^ i / + ^\ yt}l m 2i ris + j + i ) (4 . 26 ) 

00 ( Q \2n n— 1 00 / - \ fe s .\ , .\ 

B-^' r (»"+,- + 1) n<- - ') g tt r (" + k + b + 1) ("»">*) 2 v (.+«+*+y , 

where we have exploited the relation 

£ fe^=(-irnV-p), ( 4 - 27 ) 

*■ ' p=0 

which is valid for n € N + . 

We would like to make a remark at this point. It is not difficult to prove that the expression 
for A-i(s), Aq(s) and Ai(s) derived in this section reduce, in the limit of small 60, to the results 
obtained for the generalized cone [7|, [25(. In the case of the functions Ai(s) the correct limiting 
behavior for small 6q is achieved thanks to the following important property of tti(v), namely 



.(id) 



+ J2 z ^ = ° . ( 4 - 28 ) 

b=x(i) 



which can be easily proved by noticing that 

fit(l)=A(l). ( 4 - 29 ) 

This last relation comes from the fact that, for Oq — > 0, the functions fJj(f) have to coincide 
with the polynomials Di(y) appearing in the uniform asymptotic expansion of the modified Bessel 
functions. 

Let us conclude this section with a few remarks about the additional complications one would 
encounter if — A^ instead of — A^ + d 2 /4 was considered. Formally one simply has to replace 
m 2 by m 2 — d 2 /4. As a result, in the limit m — > the integrals involved do not turn into doable 
integrals from zero to infinity, but something analytically considerably more complicated evolves. 
Going back to (|3.8p it is clear that for m = 0a branch cut at k = d/2 would be present and when 
deforming the contour additional contributions due to that branch cut would arise. As stressed, 
because our focus here is the heat kernel where these compilations can easily be avoided, we do 
not discuss these additional technical complications further. 
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V. THE HEAT KERNEL COEFFICIENTS 

In this section we will use the previously obtained results for A^i(s), ^4o( s ) an d -Aj(s) in order 
to evaluate the coefficients of the heat kernel asymptotic expansion on the manifold ^#. Since the 
base manifold .jV has been left unspecified, we will find an expression for the heat kernel coefficients 
on ^# in terms of the ones on jV [7|. We will set m = in A_\{s), Aq(s) and Ai(s), because the 
mass dependence of the heat kernel is a trivial factor of e~ m . 

For the explicit evaluation of the heat kernel coefficients associated with the operator — A^ + 

d 2 /4 in arbitrary dimensions D, it is convenient to use the following formula 

[D-n\ [D-n\ , N 

^f = r (^— J ResC [-^-) . (5- 1 ) 

which is valid for n < D. By keeping the dimension arbitrary, we are able to effectively compute 
all heat kernel coefficients. In fact, the final formulas will be functions of the dimension D of the 
manifold ./#. Our task therefore is to evaluate A_i(s), ^4o( s ) an d Ai(s) about s = {D — n)/2 in 
the case m = 0. From the expressions (|4.7p . (|4.9p . (|4.19p . (|4.20p and (|4.26p in the previous section, 
we obtain 

xr(^l^)Res C ^(^^), (5.2) 

for the term Aq(s) we find 

(D-n\ , (D-n\ (sm9 ) D - n ( D - n\ ( D - n\ , , 



and, finally, for Ai(s) 

T — — Res^ — — = -(sm8 ) D - 



C 2 , l (8 ,D-n)+C 3 , l (6o,D-n)+C 4 , i (e ,D-n) 



,'D — n + i\ (D — n + C. 
< r ( J ResC^ f ^ ) • (o. I) 



Here we have defined the functions Ci(9o, D — n) as 

C(9 D n)~ l V (2fc)! V r ( J + fc + l) r ( i ^+ J ) 



^to^^U F(i + l)r(^=i) 



i=l 






C2,»(flo, g - n) = J^ Hb^- D K _ n J, - J 7 UJ j - , (5.6) 



* F ( D ~ n _|_ n\ 
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(5.7) 



i=i 



C^(e ,D-n) = ( S m9 ) n - D Y, E g ^ w p-n +i 



r(^±^ + 6 + j) 






j=l b=x(i) 
00 T-i /_. 1 I. 1 i\ n—l 

with the notation 



5», 6 = x iib (cos6» ) i+26 , Zj j6J = z iAi (cos6> ) 4+26 . (5.9) 



The formulas obtained in (|5.5p and (|5.8p can actually be simplified further. In fact, the double 
series which defines Ci(6q,D — n) can be explicitly evaluated leading to the remarkably simple 
result 

Ci(0 O) D-n)= 2 F 1 Q, ^p, ^p + 1; sin 2 O ) , (5.10) 

which is proved in Appendix[Bj Moreover, by exploiting the series definition of the hypergeometric 
function [31] 



p ( n a „.. r 2\ _ r (7) v^ r(q + m)r (/? + m) ^ 2m 

2^1 (a, £, 7, z J - — — — — - > — - — ■ - z (5.11) 

1 (q)1 (p) ^-^ to!1 (7 + m) 

and the relation (|4.27p . we can rewrite (|5.8p as follows 

x a^if-^y^^ + ^ + i + ^cos 2 ^) • (5.12) 



The relations (|5.2p through (|5.4p are all the ingredients that we need in order to write down an 
expression for the heat kernel coefficients g/n in terms of those on the base manifold. In order 
for Z{s) in (|3.19p to give no contribution to the residue of £(s) at the points s = (D — n)/2 and, 
in turn, to the value of the heat kernel coefficients, we need to set N > n — 1. In this way the 
points s = (D — n)/2 lie in the region of analyticity of Z{s). By recalling the relation (15. ip and by 
noticing that for the base JV this implies 

,jY .id—n\ n _ fd — n 



r hr Rffi C/hr . (5-13) 
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we obtain 



(sing )^ j, (sinflo)^ "^ 
20F( J D-n) Cl(0 °' jD " n) ^f r - 



("-!) 



n-1 



(sineo) D ~ n ^^(l^i) C 2 ,t(0o,I>-n)+C3,<(eo J i3-n)+C4 ) i(eo,£>-n) , (5.14) 



. 2 

2=1 

where .0^ are the coefficients of the small t asymptotic expansion of the operator — A r/ j/+(d— 1) 2 /4 

2 

on JY , and £^_i) = for n = 0. The coefficients .A^ of the heat kernel asymptotic expansion for 

2 2 

the Laplacian on the spherical suspension ^ are obtained from s^m using the relation 

If J -, /j\2fc 

fc=0 ' v 7 
The last formula, together with the explicit expression for srfn in (15.14[) . represents the main result 
of our paper. It provides a relation for the heat kernel coefficients on the spherical suspension ^# 
in terms of the heat kernel coefficients of the operator — Ajs + (d — l) 2 /4 for any smooth, possibly 
with boundary, base JV . This result is very general considering the fact that it gives any desired 
heat kernel coefficient for arbitrary dimensions D. In the leading order, the expression for A^r in 

2 

(|5.14p reduces to the one for the heat kernel asymptotic coefficients on the generalized cone once 
the limit of small 6q is considered [7]. One can obtain more explicit formulas once JV is specified. 
One might ask how the complicated 9q dependence in terms of hypergeometric functions in (|5.14p 
results from the known expressions of the heat kernel coefficients in terms of geometric tensors. This 
is in fact a direct consequence of integrating trigonometric functions raised to arbitrary powers. 
For example the leading coefficient 

Af = (47r)- D / 2 vol(^). 

For the given manifold, the volume is 

do 
vol(^) = vol(^) [ de$m D ~ 1 e = vol(^K) ^^Ci(6 ,D). 
o 

Similarly, integrating quantities involving curvature terms like for example the extrinsic curvature 



K = d cot i 



'o 



or the Riemann scalar curvature 



sin u sin 6 
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and powers thereof, results obtained from (J5.15P can be reproduced. Evaluating the heat kernel 
coefficients from the curvature tensors is actually quite cumbersome and only possible for the 



leading known coefficients 29 



30 



351 ] . For the given geometry, (15. 15ft has the advantage that any 
number of coefficients can be determined using a simple algebraic computer program. In principle 
this could be used to put restrictions on higher heat kernel coefficients, in the same way that the 
generalized cone was used [7J, |isl, |34|] . 



In the next section, we will study the situation in which the base manifold is a d-dimensional 
sphere. For this particular choice the spectral zeta function C y{s) can be written in terms of a 
linear combination of Barnes zeta functions and more specific results can be obtained [71.1121]. 

VI. A SPECIFIC BASE MANIFOLD 

In this section we will focus our attention to a particular case for which more explicit results 
for the heat kernel coefficients can be obtained . We consider the base manifold JV to be a d- 
dimensional sphere of unit radius. The eigenvalues of the modified Laplacian — A^y + [d — l) 2 /4 
are known to be 

f=(k+ d -^\ (6.1) 

where k > 0, and the eigenfunctions are hyperspherical harmonics with degeneracy 

d{k) = {2k + d-l) &+f-W . (6.2) 

The two formulas above, together with the definition of the base zeta function C/k( s ) i n (|3.4p . allow 

us to write 

-2s 



, , x v^/ , dk + d-2)\ f, d-l\~ zs , N 

M<) = !> + <*- l%3T)f { k + —) • (6 ' 3) 

3 point, it is 



fc=0 

At this point, it is convenient to introduce the Barnes zeta function, which is defined as follows 



(b(s, a\r) =} j (a + fh-f) s , (6.4) 

m=0 

valid for 9ft(s) > d, where rh and r are ci-dimensional vectors. It has been shown in [7J, [l2J that 
Cjy{s) in (|6.3p can be rewritten, after some algebraic manipulations, as a sum of two Barnes zeta 
functions 

Ck(*) = (b (2s, ^±1) + (b (2s, ^-±) , (6.5) 
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where we have used the notation £g(s,a|l) = £g(s,a). 

For the purpose of the evaluation of the heat kernel coefficients we will need the residues of 
the function C/K ( s ) at s = m/2 with m G N + . To this end it is convenient to utilize the following 



integral representation of the Barnes zeta function 



7. 



CeiB , a) = '3l^f dy ^-'\y;-\ ( 6 . 6) 

V ; 2tt J l 2^sinh d (f) V ' 

where L represents the Hankel contour. The spectral zeta function on the base manifold can now 
be represented in terms of a contour integral by using the relation (|6.5p and the formula (|6.6p . 
More explicitly, one obtains [7(] 



27r(d - 1) J L smh a i y 



CHs) = ;_7, ^ 2^- d J dy (-yf^-^r- . (6.7) 



The residues at s = m/2 for m > 2 and d > m are then easily found by an application of the 

residue theorem. They are best expressed in terms of D v defined through 

d—loo v 

= Y J Df^ V -. (6.8) 

sinh y I ^-^ v\ 

In detail one obtains 

2 m ~ d F) ~ 

ResU Vj) = (d-i)(m-md- m )r (6 - 9) 

This result can now be used in (|5,13p in order to get a formula for the heat kernel coefficients .srfiC . 

2 

One obtains 

*> _ <"-"-;> P ;-"r ( d -^) . (6.10) 

2 2 n n\{d-l)T(d-n) \ 2 J y ' 

This expression can be written in terms of a single gamma function by exploiting the formula 

„/ , n n i j « i„ (d — n + 1\ „ (d — n\ 

T(d -n) = (27T)-22 d - n -2T T — - , (6.11) 



which is a particular case of the more general product theorem for gamma functions 3l( • This last 
remark allows us to write (|6.10p as 

,^ = 2^ j d " n " 1) 5"l - (6-12) 

2 V 2 d n!(d-l)r(^f±i) V ; 

This result, applied to the general relation (|5.13p . gives an explicit formula for the heat kernel 
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coefficients £/<± in the situation in which the base is a d-dimensional sphere. In more detail we have 
(4tt)t (sing ) g -"(d-n-l) / tf-n + jA n(d _ 1} 



|5 d | ^ n!(d-l)(d-n + l) V 2 

' 2 

0F(sin# o ) D ~ n (d-n) fd-n + 2\ (d _ 1} 

n-l 
2 



2(d-l)(n-l)! V 2 



n-l 



2 ^ r«;„ a \D-n V^ ( d - n + *) /d-n + i + 2\ (d _i) 



(d-l) (sm0o) ^(n-l-i)! V 2 J^ Ul 

C2,i(eo,^-n) + C 3l i(eo,I>-n)+C4,i(eo,13-n)] , (6.13) 

where we have used the Pochhammer symbol (x) n = r(x + n)/T(x) and \S$\ represents the surface 
area of a d-dimensional sphere of unit radius 

1*1 = ^© " (6 ' 14) 

We would like to stress, one more time, that the result (J6.13J) reduces, in the leading order, to the 
one obtained in [7J once the limit of small angles 9$ is considered. The formula (16.13P allows for 
a quick computation of all the heat kernel coefficients of the manifold ^#, through the relation 
(J5.14J) . in arbitrary dimensions in the case in which the base manifold is a d-dimensional sphere. 
The resulting expressions will depend in a direct way on the functions Ci introduced in (15. 5p through 

dEHD. 

Below, we will list the heat kernel coefficients up to A.% for a d-dimensional sphere as the base 
manifold. By utilizing (15.15p . (16.13P and by introducing the function 



Ti(6 ,D -n) = C 2 , t (6 ,D - n) + C 3 ,*(#o, D-n) + C M (0 O , D-n), (6.15) 



we have 



D 



U ' r) " r Af = T^T C i(^), (6-16) 



(sme ) 1+d \S d \ ° d+1 



(4VF) " lAf = -], (6.17) 



smO.nSJ 3 4 ' 



(sJn\ A ^ = -^d^D^-^W,*-*) 



+ 4 (/ + i) sin2goC,1 ( go '- P )' (6 ' 18) 



d 



(4tt)3 . __ (d-3)(d-l) _ n Q , d 



2 



■4f = IS i " ^{Oo, D - 3) - — sin 2 O , (6.19) 



(sin# ) d - 2 |S d | § 48 zv u ' ' 16 
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D 



i47T)T -Af = (d " 5)(d "i )(5d " 3) ^o^-4) + ^ (d - 3 y- 1 ^ 1 (^D-4) 



{sme ) d - 3 \S d \ 2 1440 iV u ' ' ' v 6 

- 20F^ 3 (go,^ ~ 4) - ^ ( ^~ 3) sin^oCi^o,^ ~ 2) 

- l^sin 2 fl ^i(#o,£-2) + ^^— sin 4 fl Ci(fio,£>), (6.20) 

(4vr)i ^ (rf-5)(d-3)(rf-l)(5rf-3) (rf-3)(d-l) 

F^S = zr^z + — J- 2 {%,D-5) 



(sin0 o ) d_4 |Sd| s 5760 12 

- J*(«b, D - 5) + d2 ( d - 3 )( d ~ 1 ) sin 2 9q _ * sin 2 ^(^ D _ 3) 

iyz 4 

-^ sin 4 0o, (6.21) 

(4tt)t _^ = (d-7)(d-3)(d-l) (9 -284 + 35^) ^^ 



(sin0 o ) d_5 |'S'd| d 362880 

( d _ 5 ) (( i_3)( d _l)( 5d _3) 
- V^ ^20 Ti(9 ,D-6) 

+ V^ {d ~ 3) ^~ 1) J r 3(0 O ,D-6)-2^M0o,D-6) 

5760 
+ ^ ~ ^ ~ 1} ^BtfeoK^D- 4) -^sin 2 6 T 3 (e ,D-4) 

- d4{d ~ 3) sin 4 o Ci(0 o , I> - 2) - ^V^sin 4 0o^i(#o, D-2) 
384 lb 

J 6 
+ 384(?TI) Shl6 '''> C '"''>' D) - (6 ' 22) 

For reasons of space, the above results are expressed in terms of the functions J-i(6q, D — n) which 
are defined in equation (|6.15j) together with (|5.6p through (|5.8p . We would like to stress, once 
again, that one can easily evaluate the functions J~%(9o, D — n) explicitly with the help of a simple 
computer program. 



VII. CONCLUSIONS 

In this paper we have analyzed the spectral zeta function for a Laplace operator acting on 
scalar functions defined on the spherical suspension and evaluated the coefficients of the asymptotic 
expansion of the heat kernel for Dirichlet boundary conditions. The coefficients were given in terms 
of those on the base. More explicit result were presented once the base manifold was specified to be a 
(i-dimensional sphere. We have also provided a systematic way of dealing with the integrals defining 
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A_i(s), ^4o( s ) an d Ai(s) coming from the uniform asymptotic expansion of the Legendre functions 
used for the analytic continuation of the spectral zeta function. These represent non trivial results 
because they cannot be obtained by direct and straightforward application of previous techniques 
utilized when the radial eigenfunctions are Bessel functions. 

Since in this work we have only treated the case in which the Laplace operator is endowed 
with Dirichlet boundary conditions, the next natural step would be to extend the results to the 
more general Robin case. This analysis is possible because the uniform asymptotic expansion of 
the derivative of the Legendre functions is known 
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39l | and it has the same form as (|3.11[) . We 



expect that the analytic continuation of the spectral zeta function for Robin boundary conditions 
would follow the lines presented in this paper without any major technical complications. 

It would also be very interesting to apply the formalism developed in this paper to the evalu- 
ation of the Casimir energy for the spherical suspension following the ideas developed for conical 
manifolds in [26|, |27| . It is well known that the regularized vacuum energy is obtained from the 
value of the spectral zeta function at the point s = —1/2. One would need to analyze in detail 
the term Z{s) in (|3,19p . While for the computation of the heat kernel coefficients this term does 
not contribute to the residues of C( s )) we would be forced to deal with it in the evaluation of the 
Casimir energy. One way of expressing Z(s) would be by means of the Abel-Plana summation 
formula which has been used in [2j|. The results obtained in this way would be useful, at least, 
for a numerical estimate of the Casimir energy on the spherical suspension. 

It is not clear, at this moment, whether or not the results obtained in this work can be easily 
adapted to the evaluation of the spectral zeta function associated with — A_^. Although the 
eigenfunctions are still associated Legendre functions, this case leads to a more cumbersome uniform 
asymptotic expansion (the polynomials r2 n (^) would acquire non-polynomials terms containing the 



arcsine 33j) and new technicalities have to be considered. 
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Appendix A: The polynomials Q, n {v) 

In this appendix we list the functions VL n {y) for n up to five. In the following we will set 
7 = w/sm0Q and v is defined in (J3.12J) . By utilizing the recurrence relation (I3.15P and the 
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cumulant expansion f)3. lTj) we get 



0lM= 8-^ + (TT^)U-4 + ^) ' (AI) 



z/ 2 3i/ 4 5z/ 6 1/1 9z/ 2 3^ 4 5u 6 

1 / 1 9^ 2 9i^ 4 5i/ 6 \ ,. , 

+ tt— ^ \T- c --c7r + -o 5- , A2 



(I + 7 2 ) Vl6 16 



n , , 25z/ 3 531zv 5 221i/ 7 1105i/ 
n 3M = ^ 77777- + 



384 640 128 1152 

1 / 7 19i/ 259z/ 3 2949i/ 5 221i/ 7 1105^ 9 

+ ^[~=^ + -^r--^- + 



(1 + 7 2)3 ^ 720 32 96 640 64 1152 

1 / 7 19i/ 259i/ 3 2949^ 5 1105i/ 7 1105z/ 9 

+ (i + 7 2) 2 ^960 ~ "32" + ~~64 320 + 128 384 

1 /9^ 71i/ 3 87z^ 5 221i/ 7 { 1105i/ 9 \ 

+ (I + 7 2 ) V128 ~ ^8~ + ^6 32~ + ~W ) ' ( } 



n , , 13i/ 4 71z/ 6 531z/ 8 339i/ 10 565i/ 12 
v ; 128 32 64 32 128 

1 / 5 153i/ 2 2613z/ 4 1283&- 6 2619z/ 8 339i/ 10 565i/ 12 

+ ^7T -T7.-^^ + —^ T^T" + ^77 777" + 



(l + 7 2) 4 \^16 32 128 32 64 16 128 

1 / 5 459i/ 2 2613i/ 4 6415i/ 6 7857j/ 8 2373i/ 10 565i/ 12 
+ rr* "T77 + -77^ ttt— + —77: 777— + 



(1 + 7 2)3 ^ 16 64 64 6 4 64 32 32 

1 / 5 171zv 2 3207i/ 4 677i/ 6 2097i/ 8 3051z/° 1695z/ 12 

+ (i + 7 2)2 ^128 64" + 128 8~ + "" L6 32 + 64 

1 /9z^ 2 297^ 4 1709&- 6 3681i/ 8 1695i/ 10 565z/ 2 \ 

+ (TTt 2 ) V"64~ ~ "64" + ~ 6l 6Z - + ^32 32 J ' ( ' 

1073i/ 5 50049zv 7 186821*/ 9 44899i/ n 82825z/ 3 82825z/ 5 

5 ^ ~ 5120 7168 + 4608 512 + 1024 3072 

1 / 31 631^ 39551^ 3 539643i/ 5 3250189i/ 7 162353z/ 9 

+ (i + 7 2) 5 ^2520 ~ 128 + 768 2560 + 7168 288 

209571^ n 82825z/ 3 82825z/ 5 ' 
+ 512 512 + " 




1 / 31 1893^ 39551z/ 3 539643i^ 5 9750567z/ 7 1136471z/ 9 
+ (l + 7 2) 4 ^ 2016 + ^256 384 + 1024 7168 + 576 

209571^ n 745425i/ 13 414125z/ 15 \ 
128 + 1024 3072 J 
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1 / 31 1415i/ 65569i/ 3 116327z/ 5 10693979i/ 7 6031529i/ 9 
+ — — + 



(l + 7 2) 3 l8064 512 1024 256 7168 2304 

1302325// 1 82825i/ 13 414125i/ 15 \ 
+ 512 64 + 1536 ) 

1 /l53i/ 38503i/ 3 158319// 733859i/ 7 2476075z/ 9 972981^ n 
+ (i + 7 2) 2 I 1024 3072 + 1024 1024 + 1536 512 

579775^ 13 414125z/ 15 \ 
+ 512 1536 J 

1 fl83u 3 8613i/ 5 141923i/ 7 170509z/ 9 86067i/ n 

+ — + 



(l + 7 2 )l 512 512 1024 384 128 

248475z/ 3 414125z/ 15 \ 

512 + 3072 )' ^ ' 



Appendix B: The function Ci(9o,D — n) 

In this appendix we will provide a proof of the formula (|5.10p . Let us start with the expression 
for C\ in terms of the double infinite series (|5.5|) . 

cto 2s i- l V (2fe)i y r( J + fc + i)r( g + j-i)r(, + i) 

Cl(e ° M -A^]jr( J + i)r( s -i)r( s + i +J ^) (sm ^ ' (B1) 

n 

where we have used the fact that [31| 

In addition, one can prove by a direct calculation and by using the doubling formula for the Gamma 
function [31j 

T(2k) = 2^ mr T { ( k + ^ , (B3) 



that 

(2fc)! _ 2kT(2k) r(fc+^) 
2 2k k\ ~ 2 2k kT{k) ~ r(i) ' 

The previous results allow us to rewrite the double sum in ()Bip as 



(B4) 



By directly comparing the last expression with the definition of the hypergeometric function of two 



variables 



ml 



F, ( g , A 0, r , x. ») = r . ,^!> „, E E F( ° + '" r ffl + 1'T + - -V ■ m 

1 (all (a 1 (p) ^-^ ^-^ mlwliT + m + n 

\ / v / v / m= o n=0 ' 
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and by recalling that T(l/2) = ^/n, one finds that 

Ci(0 o , 2s) = Fi Q, s - ~, i s + 1; sin 2 O , sin 2 O ) • (B7) 

The above hypergeometric function of two variables can be written in terms of a hypergeometric 
function of one variable thanks to the relation [3l|] 

F 1 (a,/?,/? / , 7 ;x,x) = 2 F 1 (a,/3 + /3 , , 7 ;x), (B8) 

which leads to the formula (J5.10J) . once s = (D — n)/2 has been put. 
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